Non-linear transport phenomena in a two-subband system 
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We study non-linear transport phenomena in a high-mobility bilayer system with two closely 
spaced populated electronic subbands in a perpendicular magnetic field. For a moderate direct 
current excitation, we observe zero-differential-resistance states with a characteristic 1/B periodicity. 
We investigate, both experimentally and theoretically, the Hall field-induced resistance oscillations 
which modulate the high-frequency magneto-intersubband oscillations in our system if we increase 
the current. We also observe and describe the infiuence of direct current on the magnetoresistance 
in the presence of microwave irradiation. 



1. INTRODUCTION 

Nonequilibrium magnetotransport in two-dimensional 
electron systems (2DES) at large filling factors is a sub- 
ject of intense interest. Experimental and theoretical 
work includes studies of the steady-state response of 
2DESs under ac excitation in the microwave (MW) range 
and under elevated direct current (dc). In both cases, 
there appears an oscillating magnetoresistance due to 
scattering-assisted transitions of electrons between differ- 
ent Landau levels (LLs). At low temperatures, when the 
elastic scattering dominates, such transitions may occur 
due to coupling of LLs either by the MW field, owing to 
absorption and emission of the radiation quanta Hlo, or by 
the Hall field Edc which tilts the LLs. Consequentlyjthe 
microwave-induced resistance oscillations (MIROs) 
are governed by the ratio Sac = w/wc where ujc — eB/m is 
the cyclotron frequency, B is the magnetic field strength, 
and m is the effective mass of electrons. Next, the 
Hall field-induced resistance oscillations (HIROs) [l-Q 
are governed by the ratio edc = '^.Rc&Edclfu^c, where 
Rc = vf l^c is the classical cyclotron radius and vp is the 
Fermi velocity. The Hall field is defined as Edc = Ph], 
where pn = muic/e'^ns is the Hall resistivity, Ug is the 
electron density and j = Idc/w is the current density with 
Idc being the applied current and w the sample width. 

The amplitude of MIROs can increase dramatically 
with increasing MW power. In the samples with high 
electron mobility, elevated MW power transforms the 
minima of these oscillations into intervals of zero dissipa- 
tive resistance. These are referred to as zero-resistance 
states (ZRSs) The amplitude of HIROs is consid- 

erably smaller because these oscillations rely upon the 
large-angle scattering between field-tilted LLs, which re- 
quires the presence of short-range scattering potential. 
However, in the regime edc < 1, corresponding to the 
first minimum of the HIROs, the Hall field creates a 
strongly nonequilibrium electron distribution within the 
LLs, which causes a considerable decrease of the resis- 



tance This may give rise to a phenomenon of van- 
ishing differential resistance in high- mobility 2DESs 0]; 
the resulting states are called zero-differential-resistance 
states (ZdRS). The ZDRSs can also emerge from the 
Shubnikov-de Haas (SdH) oscillation maxima [l3|- The 
ZDRS phenomenon in the dc response resembles the 
ZRSs in the MW-induced response and can be explained 
with similar reasons, in terms of the negative differential 
resistivity, which leads to an instability of the homoge- 
neous current picture and spontaneous breaking of the 
sample into current domains pd| . As the instability is of 
macroscopic origin, it exists regardless of the details of 
the microscopic mechanisms of non-equilibrium magne- 
toresistance. However, these mechanisms determine the 
region of the magnetic fields where the ZRSs or ZDRS 
are observed. 

Furthermore, an interesting situation takes place when 
a 2DES is subjected to both dc and ac excitations so that 
both sources of nonlinearity are applied together. Ex- 
perimental studies (l^ . [isj have shown that the resulting 
oscillations depend on Cqc, and edc, as well as on the com- 
bination of these parameters, eac i edc- The component 
of the magnetoresistance oscillating with eac i edc can be 
described as an interference between MIROs and HIROs. 
A continuous increase of the current transforms the os- 
cillation maxima into minima (and vice versa) and may 
lead to destruction of ZRSs and to reappearance of ZRSs 
in a different region of magnetic fields. 

A unified physical picture of the nonequilibrium mag- 
netotransport phenomena described above is currently 
under development. Earlier theoretical works were fo- 
cused on the influence of MWs on the probabilities of 
elastic scattering of electrons in a system of LLs 
[T6| . which contribute to the "displacement" mechanism 
of MW-induced magnetoresistance. Later, the influence 
of both MWs and dc excitation on the electron distribu- 
tion function was found to be very important p^ . The 
deviation of this function from equilibrium is propor- 
tional to the inelastic electron-electron scattering time 
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Tin , which rapidly increases with decreasing temperature. 
Thus, the corresponding "inelastic" mechanism of non- 
linear magnetoresistance proves to be the most impor- 
tant one in the description of MW-induced magnetore- 
sistance at low temperatures, which is also confirmed 
experimentally [l8l |. A comprehensive consideration of 
possible mechanisms of MW-induced magnetoresistance 
and of their interplay at elevated radiation power and 
for different kinds of elastic disorder potentials has been 
presented in Refs. [l^ and [l^l- A theory of dc magne- 
toresistance comprising both displacement and inelastic 
mechanisms and describing the response in a wide inter- 
val of applied currents (including the region of HIROs) 
has also been developed [2l|. The interplay of both dc 
and ac excitations is studied in Refs. [l^l and [l^l, and 
the theoretical results confirm the basic properties of the 
oscillat ing magnetoresistance investigated in the experi- 
ments [l2| and [T3| . 

The nonequilibrium magnetotransport phenomena ac- 
quire qualitatively distinct features in systems with two 
populated electronic subbands, because of the possibil- 
ity of electron transitions between the two sets of Lan- 
dau levels associated with these subbands. As a result, 
even in the linear regime, such systems show magneto- 
intersubband (MIS) oscillations governed by the ratio 
A/hu!c [131 • The subband separation energy A can be 
a few meV in double quantum wells so that the MIS os- 
cillation period is comparable to the periods of MIROs 
and HIROs. Under MW excitation, the oscillating mag- 
netoresistance is determined by an interference of MIROs 
and MIS oscillations [lH, [l^l . A theory based on consid- 
eration of the inelastic mechanism of magnetoresistance 
in a two-subband system satisfactorily explains the ob- 
served oscillating resistivity as well as its temperature 
dependence in the low-temperature region [25| , , [1^ , 
[29|. Similar to single-subband systems, high-mobility 
two-subband systems demonstrate the ZRS phenomenon 
|. The non-linear response to moderate dc excitation 
30|, |3l| shows a flip (inversion) of MIS oscillations. This 
behavior signifies an inversion of the quantum contribu- 
tion to resistivity as a consequence of nonequilibrium 
electron distribution, in agreement with the theory of 
Ref. The dependence of the characteristic mag- 

netic field (where the flip occurs) on the current density 
is in a quantitative agreement with the theory (3l| . The 
response to elevated dc excitation [32[ shows the inter- 
ference of MIS oscillations with HIROs. Finally, it was 
recently found [llj that high-mobility two-subband sys- 
tems develop several regions of ZDRSs around the peaks 
of inverted MIS oscillations. 

In spite of those recent advances in studying magne- 
totransport in two-subband systems, some problems are 
still poorly investigated. In particular, neither a system- 
atic experimental study of HIROs nor their theoretical 
description has been presented. Preliminary studies of 
ZDRSs [1^ are limited to low temperatures and low dc 
currents. The behavior of two-subband systems under 
the combined action of dc and ac excitations has not been 



studied. In the present work, we intend to fill these gaps. 
We undertake studies of non-equilibrium magnetotrans- 
port in a high-mobility bilayer electron system formed by 
a wide quantum well (WQW). The high sample quality 
enables us to carry out a systematic experimental in- 
vestigation of dc transport and to describe both ZdRS 
and HIROs which appear at elevated currents. We also 
show and discuss the influence of a dc electric field on 
MW-induced ZRSs at a chosen frequency of 143 GHz. A 
supporting theoretical consideration, based upon a gen- 
eralization of existing theories to the two-subband case, 
reasonably explains the different regimes of non-linear 
magnetotransport . 

The paper is organized as follows. Section II briefly 
presents experimental details. In Sec. HI we describe 
the magnetoresistance measured in a wide range of ap- 
plied dc excitation and compare our data in the regime 
of HIROs with the theoretical calculations. In Sec. IV 
we analyze the influence of moderate and strong dc exci- 
tations on the magnetoresistance under MW irradiation. 
Concluding remarks are given in the last section. 



2. EXPERIMENTAL SETUP 

For our experiments, we use high-quality WQW's with 
a well width of 45 nm, high electron density ~ 
9.2 X 10" cm-2, and a mobility of // ~ 1.9 x 10*^ cm^/V s 
after a brief illumination with a red light-emitting diode. 
Several samples in Hall bar geometry (length / x width 
If) of (500 X 200) /xm^ have been studied. The electrons 
in our system occupy the two lowest subbands with a 
subband separation of A = 1.4 meV, extracted from the 
periodicity of low-field MIS oscillations. Measurements 
have been performed between 1.4 and 4.2 K in a cryo- 
stat with a variable-temperature insert. For experiments 
under MW irradiation, a waveguide is employed to de- 
liver linearly polarized MW radiation down to the sam- 
ple. We record longitudinal resistance using a current of 
1 iiA at a frequency of 13 Hz in the linear regime. Di- 
rect current Idc was applied simultaneously through the 
same current leads to measure the differential resistance 
r^x = dVxx/dldc- In the next section we describe the 
measurements of the resistance under dc excitation and 
its theoretical analysis. 



3. MAGNETORESISTANCE UNDER DC 
EXCITATION 

In Fig. [l]we present the differential resistance Txx as a 
fimction of the magnetic field for Idc=Q, Idc~50 fJ-^, and 
l]jc~^^0 For ldc=0, the magnetoresistance exhibits 
well-developed MIS oscillations and thereby confirms the 
existence of two populated subbands It should be 

noticed that the ac current of 1 ^A is already high enough 
to invert MIS osciUations for B < 0.2 T. At 1.4 K, we 
also observe SdH oscillations which are superimposed on 
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FIG. 1: (Color online) Differential resistance r^^ as a func- 
tion of the magnetic field at 1.4 K. ZdRS are present for 
high DC current, see arrows. In addition, HIRO's appear 
for /dc=150 /iA (see the box) at low magnetic fields. Inset: 
measurement setup for differential resistance. 



the MIS oscillations for B > 0.25 T. We first apply a 
moderate current of 50 /iA and find that the MIS oscilla- 
tions are inverted in the whole range of magnetic fields. 
Around 0.27 and 0.4 T, the inverted MIS peaks evolve 
into wide regions of ZDRSs. We also find a splitting of 
the oscillation maxima, e.g., at 0.23 and 0.32 T similar 
to that in dc-driven non-linear magnetotransport in dou- 
ble quantum wells using a high current (3l| . A further 
increase in the dc electric field, see trace for Ij^c=150 fiA, 
leads to vanishing ZDRSs at B=Q.27 T and to disappear- 
ance of the splitting. Apart from that, the most intrigu- 
ing features appear around _B=0.1 T, where the inverted 
MIS peaks are transformed into enhanced ones. Other 
MIS peaks, e.g. at 0.2 T, remain inverted. We iden- 
tify the feature around 0.1 T with the first HIRO peak. 
Indeed, if Jdc= 150 fxA the condition edc = ^, which corre- 
sponds to the maximum of the first HIRO peak, occurs at 
B ~ 0.07 T, not far from the point B ~ 0.09 T where we 
actually observe the maximum. The shift of the observed 
maximum to higher fields occurs because the ocsillating 
contribution to resistivity is exponentially suppressed by 
the disorder in the region of low magnetic fields. 

We now focus on HIRO's at elevated currents. Figures 
[2] (a) and [2] (b) give an overview of magnetoresistance 
under dc excitation. In this range of magnetic field from 
0.05 to 0.15 T and at a temperature of 1.4 K, MIS oscil- 
lations are strongly infiuenced by a dc electric field. The 
SdH oscillations, which occur at higher magnetic fields 
are damped due to electron heating. For Idc between 150 
and 200 /xA, we first see that all MIS oscillations that 
were inverted for 50 /iA, are now strongly enhanced at 
B > 0.07 T. In addition, the features around 0.06 T in 
Fig. [2Ja) become inverted once again with increasing dc 



excitation. A further increase in the dc excitation from 
225 to 300 /iA modifies the differential magnetoresistance 
in a similar way, by modulating the MIS oscillations ac- 
cording to the HIRO periodicity. This behavior is also 
illustrated in Fig. [2jc), where we plot r^x for both neg- 
ative and positive B. In order to study MIS oscillation 
behavior with increasing electric field, we mark three MIS 
peaks (I) to (HI), which we investigate in detail in Fig. 
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FIG. 2: (Color online) (a), (b) Normalized differential resis- 
tance rxx{B) /rxx{Q) as a function of the magnetic field for 
several chosen dc currents showing HIRO's for elevated Idc- 
(c) Differential resistance r^x for 150 and 300 /iA. Several MIS 
oscillations (fixed B) are marked with (blue) lines and named 
from (I) to (III). No SdH oscillations are present in the chosen 
range of magnetic fields. 

In Fig. [31[a) we first show r^x as a function of Idc 
at 5=0.269 T (as well as the corresponding voltage 
Vxx) where we observe a ZDRS. This state is created 
at /^™ ~40 /tA and is observed until a maximal current 
jmaa: ^^^^ Figure[3l^b) presents differential resistance 
as a function of Idc for MIS peaks (II) and (HI) with the 
corresponding magnetic fields indicated. For 0.07 T (HI), 
the resistivity shows oscillations which are clearly peri- 
odic with the current, and the period corresponds to that 
of HIROs. 

In Fig. |3]we show the normalized differential resistance 
as a function of i? at a temperature of 4.2 K. The be- 
havior is basically the same as in Fig. [21 For Ioc=^i we 
observe MIS oscillations whereas SdH oscillations are not 
present anymore. The amplitudes of HIRO's are weak- 
ened at higher temperatures due to enhancement of in- 
elastic electron-electron scattering and the corresponding 
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FIG. 3: (Color online) Differential resistance as a function of 
the applied DC current Idc for the inverted MIS peak which 
develops into a ZdRS (a) and for two chosen MIS peaks II 
(-8=0.096 T) and III (5=0.07 T from -400 to 400 /iA) at 
lower magnetic fields (b). 



decrease of the quantum lifetime of electrons [20|, |3j 




FIG. 4: (Color online) HIRO's at a temperature of 4.2 K: 
a) for /dc=0, 50, 100, and 150 /lA and (b) for 200, 250 and 
300 /^A. The oscillating response to dc excitation weakens 
with increasing temperature. SdH oscillations are not present 
in the chosen range of magnetic field at this temperature. 

Let us analyze the observed behavior in more detail. A 
theory which satisfactorily describes the ori gin of ZDRSs 
in our samples has been presented in Ref. [33|. For the 
regions of magnetic field around 0.27 and 0.4 T, corre- 



sponding to inverted MIS peaks, the differential resis- 
tivity calculated within the assumption of homogeneous 
current appears to be negative if T is low enough. The 
negative differential resistivity, according to Ref. [Tlj . 
leads to current domains and to ZDRSs in the observed 
magnetoresistance as in the case of single-subband sys- 
tems @, . The theory [s^ also explains the splitting of 
the maxima between inverted MIS peaks as a result of the 
influence of Landau quantization on inelastic electron- 
electron scattering in two-subband systems (see also Ref. 
[lH). The effects of higher currents described above are 
consistent with this theory. In particular, the ZDRS is 
suppressed owing to heating of the 2DES and the related 
decrease of the quantum lifetime of electrons, which leads 
to a suppression of all kinds of magnetoresistance oscilla- 
tions, especially for lower magnetic fields. The splitting 
disappears because at elevated currents (e^c ^ 1) the 
infiuence of electron-electron scattering on the distribu- 
tion function of electrons becomes inessential [2l[, this 
function is stabilized by the Hall field rather than by the 
inelastic scattering. We have confirmed these properties 
by a direct numerical calculation of the magnetoresis- 
tance. By including the short-range scattering potential 
into the theory (see the details below), we also checked 
the presence of the first HIRO peak near S = 0.1 T at 
Idc — 150 ^P^■ This peak coexists with the ZDRS around 
0.4 T, in full accordance with our experimental data (see 
Fig. [1]). 

We now concentrate on theoretical calculations of mag- 
netoresistance under strong dc excitation and weak mag- 
netic fields, in the regime where HIROs are observed. 
In the description of electron transport, we assume elas- 
tic scattering of electrons by a random static potential 
of impurities or other inhomogeneities. If the region of 
magnetic field below 0.3 T we apply the approximation 
of overlapping Landau levels. The quantum contribution 
to resistivity in this regime is proportional to the square 
of the Dingle factor d = exp(— Tr/rqWc), where Tq is the 
quantum lifetime of electrons. Under these conditions, 
non-linear magnetoresistance in single-subband systems 



is described by the theory of Ref. [21[ . A generalization 
of this theory to two-subband systems is straightforward 
since the subband separation A is small compared to the 
Fermi energy £i^, so the difference in Fermi velocities and 
scattering rates for different subbands can be neglected. 
In the regime of classically strong magnetic fields, the 
resistivity of the electron system is found from the fol- 
lowing relation between the density of the applied current 
j = /dc/'w^ and the longitudinal electric field i?| I = Vxx/l'- 



jpo|l - rf-^ 1^1 + cos ■ 
2(97(0/50' 



27rA 



90 Ttr/Ti„+Ttr/Tq-j{Q 



(1) 



where po — m/e'^nsTtj. is the classical Drude resistivity, 
Tf j. is the transport time and r^n is the inelastic scattering 
time. The dimensionless function 7(C) is given by the 
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expression 
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where Jfe(x) are the Besscl functions and 1/t((?) is the 
angular-dependent scattering rate at the Fermi energy. 
The integral in Eq. ([2]) describes the averaging over the 
scattering angle 9. The second part of Eq. © shows an 
equivalent representation of 7(C) in terms of angular har- 
monics of this scattering rate [21|. The transport time 
Ttr is defined in the usual way, as the inverse of the angu- 
lar averaging of the factor (1 — cos 9) /t{9). The dimen- 
sionless parameter C — nedc characterizes the Hall-field- 
induced tilt of Landau levels. One can write C through 
the current density as 



c 



/ 47r3j2 



(3) 



The formal difference of the theory described by Eqs. 
^ from the single-subband theory of Ref. [2l| is given 
by the MIS oscillation factor [1 -|-cos(27rA/?ia;c)]/2 in the 
quantum contribution to the resistivity and by an extra 
factorof l/V2inEq. (3). 
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FIG. 5: (Color online) Normalized differential resistance at 
T = 4.2 K for Idc = 300 and Id„ = 200 /xA. Broad lines 
(black) are the experimental plots while narrow lines show the 
calculated magnetoresistance for the cases when the relative 
contribution of short-range scatterers into the transport scat- 
tering rate l/rtr is 50 % and 90 %. The amplitudes of the 
oscillations increase with increasing content of the short-range 
scatterers. 

The calculation of the differential resistance from Eqs . 
dJl-dS]) is based upon the model of mixed disorder |20l - |22l |. 
We assume that the random static potential contains two 
components: the short-range and the long-range ones. 
The corresponding scattering rate is expressed as 
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where the angle-independent Ts gives the contribution of 
the short-range potential into the total scattering rate. 
The long-range contribution is approximated by an expo- 
nential form Tj~^ (X exp(— Zc?), where hq = 2ppsin[9 /2) 
is the transferred momentum, pp is the Fermi momen- 
tum, and Ic is the correlation length of the potential. For 
numerical calculations, we have used Ttr found from the 
zero-field mobility and determined the quantum lifetime 
at r = 4.2 K as Tq = 6 ps, based on the comparison 
between the measured and calculated magnetoresistance 
in the linear transport regime. The inelastic scattering 
time is estimated [131 '''in — ^£f/T'^, where Tg is the 
electron temperature. The appreciable heating of the 
electron gas by the applied field is obtained from the en- 
ergy balance equation taking into account the energy loss 
due to electron-phonon interaction (for example, the cur- 
rent of 300 /iA is found to heat electrons from T = 4.2 
K to Te ~ 4.8 K). In Fig. 5 we plot the calculated nor- 
malized magnetoresistance for two different values of the 
applied current Idc and compare it with the experimen- 
tal magnetoresistance. While for both currents we have a 
good agreement in the HIRO periodicity, the amplitudes 
of the oscillations expected from the theory are smaller 
than the experimental ones in the low-field part of the 
plot. Naturally, since the HIROs exist due to backscat- 
tering processes, their amplitudes increase with increas- 
ing content of the short-range scatterers. This also leads 
to a better agreement between the theory and the exper- 
iment. Therefore, based on the large HIRO amplitudes 
observed, we may conclude that the amount of short- 
range scatterers in our samples is significant. 



4. INFLUENCE OF DC EXCITATION ON THE 
MW-INDUCED MAGNETORESISTANCE 

In this section, we present experimental results of com- 
bined action of dc and ac excitations in our bilayer sys- 
tem. As it has been shown in Ref. [l^li a ZRS occurs 
around B ~ 0.27 T for 143 GHz at a temperature of 
1.4 K. We focus on that particular frequency with the 
corresponding MW electric field of ~ 4.2 V/cm for 
dB attenuation. As in previous studies [H, [23, [28| . 
MIROs and ZRSs are not sensitive to the orientation of 
linear polarization. The immunity of MIROs and ZRSs 
to the sense of circular polarization has also been found 
in Ref. [ssj where the authors employed a quasioptical 
setup in which linear and any circular polarization can 
be produced in situ. 

In Fig. ini we plot the differential magnetoresistance 
for /(jc=0 and for different DC excitations. The zero-bias 
[Idc = 0) trace shows a ZRS developed from the inverted 
MIS peak at 5=0.27 T [dashed-dotted line (1)], the other 
inverted MIS peaks at 0.15 [line (3)] and 0.1 T, as well as 
strongly enhanced MIS peaks around 0.19 T [line (2)] and 
0.4 T. If we apply a dc bias starting with /rfc=40 /xA, we 
obtain a strongly modified oscillation pattern: the MW- 
enhanced MIS peaks exhibit strong damping or even in- 
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FIG. 6: (Color online) Magnetoresistance under an AC exci- 
tation of 143 GHz {Idc~0) and under combined ac and dc for 
several chosen DC currents of /dc=40, 150 and /oc=220 /iA. 
The ZRS at 0.27 T is destroyed for 150 and 220 fiA. A ZDRS 
appears under MW excitation from the enhanced MIS oscil- 
lation at 5=0.4 T for Idc >40 fiA. Inset: r^^: under MW 
excitation for fixed B for MIS peaks from (1) to (3) as a func- 
tion of Idc (dashed line without MW's, see also Fig. [3ja)). 



version, while the effect of Idc on the MW-inverted MIS 
peaks (including the one transformed into ZRS) is very 
weak. A further increase in dc excitation has two conse- 
quences: (i) the ZRS around B = 0.27 T vanishes and 
(ii) a ZDRS develops under MW irradiation appears in 
a wide region around B = 0.4 T. For the highest DC 
excitation of Idc = 220 /iA, we also find well-developed 
HIRO's for B < 0.2 T. The dependence of the magne- 
toresistance on the dc current is non-monotonic for the 
traces (2) and (3). We have also carried out B sweeps 
at different Idc up to 220 fiA (not shown here) and mea- 
surements of the resistance at constant magnetic field as 
a function of Idc (see inset of Fig. For the ZRS at 
5=0.27 T, we find that a dc excitation leads to a van- 
ishing ZRS for Idc > 50 /iA which corresponds to an 
electric field Edc > 37 V/m. At B = 0.41 T, a ZDRSs 
develops under dc excitation develops for relatively low 
dc currents, Idc — 28 fiA giving rise to Edc —39 V/m. 

A remarkable fact is that we sec the ZDRSs in the 
same region of fields, B ~ 0.27 T, regardless of the pres- 
ence or absence of ac excitation, although with AC ex- 
citation we see the ZDRSs even in the linear response. 
As seen from the inset to Fig. [SJ this state vanishes for 
Idc > 50 /iA. For comparison, in the case of pure dc exci- 
tation the ZDRS appears at Idc > 40 /xA and vanishes at 
Idc > 80 fiA, see also Fig. EI a). The observed behavior is 
imposed by the ability of both kinds of excitation to cre- 
ate in this region of magnetic fields a partially inverted 



electron energy distribution which leads to negative re- 
sistance. Thus, our observation is consistent with the 
domain model for the origin of both ZRSs and ZDRSs. 

The basic features of the observed behavior can be un- 
derstood in terms of interplay of the effects of ac and 
dc excitations. The MW excitation cither enhances or 
inverts MIS peaks (or groups of such peaks) in accor- 
dance with MIRO periodicity [25|: the enhancement oc- 
curs around Cac = n + 3/4 (MIRO maxima) while the 
inversion corresponds to eac = n + 1/4 (MIRO minima) , 
n being an integer. The DC excitation tends to invert all 
of the MIS peaks in the region of higher magnetic fields, 
where e^c < 1- Thus, the MW-enhanced MIS peaks 
should be suppressed and eventually inverted by the DC 
excitation, while the peaks which arc already inverted 
by MWs arc expected to be little affected. However, 
the rapid suppression and inversion of the MW-enhanced 
MIS peaks already at a moderate DC excitation suggests 
that the interplay of ac and dc excitations is not reduced 
to a pure superposition of their effects. To make this 
point clear, we apply a simple theory valid at ^ 1 
and low temperatures, when the inelastic mechanism of 
non-linear response dominates, and use the approxima- 
tion of overlapping LLs. The non-linear conductivity of 
a single-subband system in this case is given by Eq. (15) 
of Ref . [l3l and can be straightforwardly generalized for 
a system with two closely spaced subbands. Hence, the 
differential magnetoresistance in written in the form 
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cos- — F{eac-,edc), (5) 
nujc J 



F 



(1 + Df - 10(1 + D)Q - 3Q^ + A{Q-l-D) 
[l + D + Qf 



where A = 2neac'P^ sm2'Keac, D = V^sin^neac, Vlu oc 
TinE^ is a dimcnsionless quantity proportional to MW 
power [l3|, and Q — (Ti„/2Ttr)(7redc)^ is a dimcnsionless 
quantity proportional to /^^. In the MIRO minima and 
maxima, A = ±27reac^tjj respectively, while D = Vui/'i, 
in both cases (notice that |^| ::g> D for the relevant region 
bj > LOc). For negative A (MIRO maxima, enhanced MIS 
peaks), when all the Q-depcndent terms in the numer- 
ator of F have the same sign, the decrease of rxx with 
increasing Idc goes much faster than in the absence of 
MWs. Physically, this is related to a stronger effect of 
the dc field on the electron distribution if the latter is 
already modified by MWs. This explains the dramatic 
suppression of the magnetoresistance peaks by the dc ex- 
citation. For positive A (MIRO minima, inverted MIS 
peaks), the MW term (cx A) in F considerably slows 
down the effect of Idc on magnetoresistance: the electron 
distribution is modified by MWs in such a way that it 
is only slightl affected by the dc field. However, in both 
cases the increase in Idc leads to an additional heating of 
the 2DES. As a result, both and decrease, and so 
does the amplitude of oscillations. In Fig. [6] we see this 
effect not only in the vanishing ZRS at 0.27 T [line (1)] 
but also in decreasing of the amplitude of the inverted 
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peak at 0.19 T [line (2)]. We have estimated that the 
current /^jc = 150 /^A leads to an additional heating of 
the 2DES in our sample by approximately 1 K, which 
may explain the observed behavior. Below 0.2 T, Eq. 
(5) is no longer valid for the highest Idc used. In this re- 
gion of fields, the non-linear magnetoresistance can not 
be simply described by the ac and dc response of the 2D 
system (superposition) but rather as an interplay (inter- 
ference) of MIROs and HIROs. 

5. CONCLUSION 

We have studied magnetotransport in a high-quality 
bilayer electron system under a strong dc excitation. We 
have found wcll-dcvcloped HIROs which modulate the 
MIS oscillations in our bilaycrs. A theoretical calculation 
of the non-linear magnetoresistance gives a good agree- 
ment with experiment concerning the shape and period- 
icity of the oscillation picture. The comparison of theory 
and experiment indicates a significant role of electron 
scattering by the short-range random potential in our 



samples. Furthermore, we show that HIROs in the re- 
gion of low magnetic fields coexist with zero differential 
resistance at higher fields. 

We have also examined non-linear transport under 
both ac (microwave) and dc excitations. Even a mod- 
erate dc excitation strongly modifies the MW-induced 
resistance. This effect is consistent with a theoretical 
model predicting that the combined action of two kinds 
of excitations is not merely a superposition of ac and dc 
effects. It is found that for the chosen MW frequency 
both these excitations lead to ZDRSs in the same region 
of fields, around B = 0.27 T. The properties of ZDRSs 
with and without MW excitation are in accordance with 
the domain model which explains the vanishing resistance 
as a result of instability of homogeneous current flow un- 
der conditions of negative differential resistivity. A better 
understanding of the interplay of ac and dc excitations is 
desirable and requires further experimental and theoret- 
ical work. 
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(Project number Ph 109/08), FAPESP and CNPq 
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